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Abstract 



Critical points of semiclassical expansions of solutions to the dispersionful Toda hierarchy are 

C ' considered and a double scaling limit method of regularization is formulated. The analogues of the 

critical points characterized by the strong conditions in the Hermitian matrix model are analyzed 

and the property of doubling of equations is proved. A wide family of sets of critical points is 

introduced and the corresponding double scaling limit expansions are discussed. 
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1 Introduction 

In a recent work pLj we considered the large N expansion of the Hermitian matrix model 

^ ^ ;, \ 1 



(1) 



k>\ 



from the point of view of the dispersionful Toda hierarchy [2] 

e— = [(£^)+,£], C = A + u + vA~\ A := exp (e d,) , e:=g/N, 



(2) 



where ()+ denotes the polynomial part in A. If the equilibrium measure in the large A^ limit is 
supported on a single interval (one-cut case), then the partition function becomes a tau- function 
T{x,t) of this hierarchy at x = g. The corresponding solution of ^ is characterized by two string 
equations which can be solved in terms of a semiclassical expansion 



J^e'^nWlx.t), t; = ^62^'7;(2'=)(x,t) 



(3) 



fc>0 



fe>0 



The coefficients {u^ \v^ ') can be expressed as rational functions of the leading terms {u^^> ,v^^'] 
and their derivatives with respect to x. 

A method to determine the expansions ([3]) was provided in [l]. It uses two functions 



E Rk{u,v) J ._sr' Tkiu,v) 



i?( 



— -to 



1, 



k>0 fc>0 

which are related to the resolvent TZ := (z — C)~^ oi the Lax operator C through the equation 

21 



(^-IT 



A 7^, 



(1+ T) 

so that they generate the trace densities Res(7?-A^-') for j > 0, where Res(^Cfc A*') := cq. These 
functions are determined by the equations 



^T + T[_i] + -(u-z): 



0, 



(4) 



T2 



41] 



1. 



where we are denoting 



F[^]{x) :=F{x + re). 
As it was shown in [Ij the string equations for (u, v) reduce to the system 

dz 



2Tri z 



dz 



V,iz,t)Riz,u,v) =0, 



(5) 



-—Vz{z,t)T[_i]{z,u,v) 



KJ-y 



27ri 



-2x, 



where V{z,t) := YlTLi ^j -^"'i ^^^ 7 is a large positively oriented closed path. For e = the system 
dS]) reduces to a pair of hodograph type equations 

/.(t,n(°),r;(0)) = 0, /4t,nW,«W) = 2x, (6) 

where 

f{t,u,v) := (b dzV^{z,t) ({z - uf - Av)^ = -2vti - 4uvt2 + . . . . (7) 

These equations determine {u^^',v^^') and represent the planar limit contribution to the partition 
function of the Hermitian model [6j. However, near critical points of ([6]) the functions (u^^' ,v^^') 
are multi- valued and have singular x-derivatives {gradient catastrophe), so that the semiclassical 
expansion ([3]) breaks down. In terms of the Hermitian matrix model this situation corresponds to 
the critical points of the large N expansion of .^Ar(t) and, as it is well-known [3]- [12], a double scaling 
limit method of regularization is available. This method leads to an important nonperturbative 
approach to two-dimensional quantum gravity. 

Recent research |13)-|16) shows that the double scaling limit expansions are also relevant in the 
asymptotic analysis of solutions of dispersionful integrable systems. They supply good approxi- 
mations in the transition regions from the semiclassical to the oscillatory regimes after the time of 
gradient catastrophe. These applications motivate the interest in reconsidering the double scaling 
limit method from the point of view of the theory of string equations for integrable systems of dis- 
persionful type. The formulation of the double scaling limit within the theory of the dispersionful 
Toda hierarchy was already addressed in |lj, but only the case corresponding to u = and even 
potentials V{—z) = V{z) was analyzed. The purpose of the present paper is to analyze the general 
case that exhibits a much richer structure of critical behaviours (see |10j-|12) for its relevance in 
the Hermitian matrix model of random surfaces and quantum gravity). 

According to the strategy of the double scaling method we look at the systems dl])-® as 
singular perturbation problems in the small parameter e and consider the semiclassical expansions 
near critical points as outer expansions of the solutions of ([4])-([5]). Consequently the problem 
consists in characterizing appropriate inner expansions 

x = xc + e2"x, u = nc + ^eV^\ v = Vc + ^e^v'^^\ e = e2"+^ (8) 

i>2 j>2 

with coefficients depending on stretched variables {x, t) . This type of expansions were determined in 
the theory of random matrix models through several methods of analysis of discrete string equations 
for orthogonal polynomials in the large A^ limit. The present work provides an alternative scheme 
based on the characterization of solutions 



E 



€^R^^\ T = Y,^^T^^\ 



of the resolvent trace equations (j3]) of the dispersionful Toda hierarchy. This scheme establishes 
the existence of the series ([8]) to all orders in e and leads to the characterization of the sublead- 
ing coefficients {w'^' ,v^'^') in terms of ordinary differential equations involving the Gel'fand-Dikii 
polynomials of the resolvent trace expansion of the Schrodinger operator. 

We analyze two different classes of critical points. Firstly, we consider the sets C„ of critical 
points which underlie the Hermitian matrix model of quantum gravity [TQ]-[11] and we determine 



inner expansions ^ which satisfy the doubling property 

with Gn,± being the Gel'fand-Dikii polynomials. Secondly, we characterize a wide family of sets 
Cimn of critical points where I > 1, m = 1, . . . ,1 and n > 2, describe the branching behavior of 
{u^^\ v^'^') near the critical points. Then we develop the double scaling limit method for the case 
n = 2 with stretched variables 

X = Xc + e^x, tj = tcj + (^ij, j > I, e = e^'^. 

For Cii2 we find that the subleading coefficients satisfy 

with v^'^' being determined by the Painleve I equation. For the other cases Cim2, / > 2, we find 
again the doubling property. Indeed, the functions 

turn out to be determined by a pair of decoupled Painleve I equations. 

The double scaling limit expansions at critical points of the types Cimn can be applied to the 
asymptotic analysis of solutions of the dispersionful Toda hierarchy. As an illustrative example 
we end our work by applying our results to regularize some critical processes of ideal models of 
Hele-Shaw flows P^j-[20] based on the Toda hierarchy. 

Our scheme can be formulated with other scaling ansatze for stretched variables like those of 
Galilean type used in [13]-[16]. Furthermore, it can be applied to the study of the double scaling 
limit of the normal matrix model which is related to another solution of the two-dimensional 
dispersionful Toda hierarchy [2T]-[23] determined by a pair of string equations. We also expect 
the methods of this paper to be applicable to dispersionful versions of the multi-component KP 
hierarchies which arise in random matrix theory and Dyson Brownian motions |26j . 

The layout of the paper is as follows: in section 2 we show how inner expansions of the gen- 
erating functions M and T can be determined from the resolvent trace equations ^. In section 
3 we characterize the double scaling limit expansions at critical points satisfying strong condi- 
tions. Finally, we devote section 4 to introduce families Cimn of critical points and to develop the 
corresponding double scaling limit method. 

2 Inner expansions of the resolvent functions 

Let us consider the hodograph system ([6]) 

/, = 0, fu = 2x. (9) 

Notice that the function / defined in ^ satisfies 

fuu -V fvv = 0. (10) 



The hodograph system determines a solution of the dispersionless Toda hierarchy provided that the 
imphcit function theorem apphes for obtaining (u{x,t),v{x,t)) from ([9]). For example ([9]) implies 

Uti = Vx, Vti = VUx, 

so that 

Vxx - {logv)t^ti = 0, 

which is the hyperbolic version of the long wave limit of the Toda lattice equation. It is also easy 
to see that the solutions of Q verify 

which is the dispersionless limit of the focusing nonlinear Schrodinger equation [16j . 
The set of critical points {xc,tc,Uc,Vc) of ([9j) is characterized by the system 

Our main aim in this work is to characterize inner expansions for solutions {u, v) of the system 
(j4|)-([5]) near critical points of ^. 

2.1 Symmetric variables 

Let us start by analyzing the resolvent trace system ([4]). It is form invariant under the change of 
variables 

x = xe + e2"x, e = e2"+\ n > 2. (12) 

Indeed since edx = e c^ we can rewrite dH as 



T + Tpj^ + -{u - z)R = 0, 



t2_^^^mm^ = i, 



(13) 



where we are denoting 

F[^](x) ■.= F{x + re). 

We will prove that if {u, v) are expansions of the form 



u = Uc + '^ e^u^^\ V = Vc + ^ e^v^^\ (14) 

j>2 j>2 

then the solution of (J13p can be expressed as 

T = ^ p¥^^ , M = ^ PR^^\ (15) 

i>o j>o 

where the coefficients {T , R ) are differential polynomials in the coefficients of the expansions 
of {u,v). 



An important property of the combined system dl])-© is its invariance under the transformation 

u{e,x) :=n(— e, x + e), v{e,x) := v{—e,x), 



T{e,z,x) ■.= T{-e,z,x + 2e), M(e,z,x) := M(-e, z,x + e). 
Thus it is convenient to introduce the symmetric variables 

w{e) := u{€,x- -) = upY/2]i'^)' 

U(e) :=T(e,z,x-e) =Tp^(e) Y{e,z,x) ■.= R{e,z,x - ^) = R^-^^^{e), 
which transform under ()16p in the simple form 

w{e) = w{-e), U(e) := U(-e), V(e) = V(-e). 
In terms of {w,v,l],Y) the system (fT3|) reads 

'U[^+U+;(^[I72]-^)V[T72]=0, 

Moreover, from ()20p one easily deduces the following linear equation for U 

{z - w^-^Yj2]) (^ - ^[172]) (^ - ^[372]) (U[T] - U) = V[T] (^ - -"^[^172]) (U[2] + U[^) 



T;(z-u;[372]))(Upjj+U). 



It can be expressed as 



(5^e'=Jfc)u = 0, 



fc>i 



(16) 



(17) 
(18) 

(19) 



(20) 



(21) 



(22) 



where J^k are linear differential operators depending on the coefficients of the expansions of (w, v) 
of the form 

Jk = {z - Uc) Vc Jkfl + {z- Uc) Jk,l + re Jk,2 + Jk,3, 

where 

re = reiz) := {z - Ue)"^ - 4Ve, 

and J7fc i are z-independent. For example, the first few are 



Jl = {Z- Uc) re ftr , J2 = 7;iz- Ue) re 9|, 



(23) 



Ji 



Z -Ue 



a| - 3rcn(2) d^-{z- Ue) {vedir + 4t;(2) cW + 24^0 - "^Ve (2u(2) a^ + J?)^_ 



' c^x 



We will now use both (|2U|) and (|21|) to determine the expansions of U and V. 



If we substitute in ([20]) the expansions 

w = Uc + y^e^w^-'' , V = Vc + y ^e'w-'', 
i>2 i>2 

j>0 j>0 

then by identifying the coefficients of the powers of e we find 

Z!!-J_ 77(o)_i^^ Z^_TT«-n 

7 ~ 1/2 ' ^ 1/2 ' 7 ~ ' 

and the recurrence system (j > 2) 

fc+i=j, fc>l k+i=i,k>l ^ ' 



*1 + «2 + fcl + ^2 = j 
ii > 2 



„ V^'^^ ' 



'^ 7. 7. 9 Z^ 



z z 2 

i+k=j, l<*,fc<i— 1 



n + «2 + fci + fe = i 

1 < ii , 22 < J — 1, A,'! + ^2 even 









+ 






n + 22 + fci + A;2 + i = J 
Z > 2, fci + A;2 even 

-(i) v^'^ 

These equations are Hnear with respect to U and . Furthermore the determinant of their 

z 

corresponding coefficients is Tc . 



Thus one obtains 



1 t7(2) 1 

17(2) = J^[4t;,iu(2) + 2{z - u,)v^% ^— = ^[2U(2) + (z - n,)wj(2)], 

2 ■S^ 2 



77(3) _ 1 ,,„. , 



(3) 



J/^-^' = -:^[4vJij('V + 2(z - Uc)U 



y{3) -I 

77(3)1 J^ — ^ — [2I;(^^ + (^ - ne)u;(^)], 



z 2 



3 
2 

re 



4„^^(4) + l^g) + 4^(2)^(2) + 2(z - n,)u(^) 



+ 



1 

+— 



¥'^ 1 



-2 2 



2^(4) + l^g)+^(2)^ + (^_„^)^Zj(4) 



+ 



re 



2vjj^^ + eW^^' + 6t;(2)' + {z- Uc) (vjv'^ + eu^^^uJ^^)) 



and by using induction one proves that the coefficients of the expansions of U and V are of the 
form 

7j(2i) ^ V- aj '{z-Uc) + Pl ^ Tj(2i+i) ^ V- aj ' {z - Uc) + Pj ^ ^ 



1=1 






1=1 



J+-2 



(24) 






1=1 



re 



1=1 



'+5 
re 



Moreover by substituting (|24p in the hnear equation (|2ip and taking into account (j23p it can 
be proved by induction that the fi 

(u(2),...,I;(^-2'+2),lZj(2)^...^^(.-2Z+2))^ 



be proved by induction that the functions aj , /3|* , 7^ and rji are differential polynomials in 



2.2 Invariant solutions 

The solution (U,V) of the system ()20p is uniquely determined by (tfjf). Hence if we assume 
that {w,v) are even functions of e, then as a consequence of the invariance of (j20p under the 
transformation (116p we deduce that (U, V) are even functions of e too. In this way we have 



+ y2f'2 3-i^f) 



w = Uc+ y ^ w 
J>1 



V = Vc+ > i^^V^^^^ 



j>0 j>0 

These solutions will be henceforth called invariant solutions of the resolvent trace equations. 



3>l 



(2j) (2i) 

As we have seen in the above subsection the coefficients U and V can be expanded in 
powers of re with leading terms 

(y ., ^^(2j) o(2i) -j^(2j) / „ \^,(2i) (2i) 

^'"- ' ''^, '' ^°(^)- V- ' ''l/^ -"{-h\ <-) 

Furthermore, by identifying the coefficient of e^^ in the first equation of the system ([20]) and by 
taking into account that {z — u^)'^ = Tc + 4fc it follows that 

Now we prove the connection between these coefficients and the Gel'fand-Dikii differential polyno- 
mials of the KdV theory. 

Theorem 1. The functions 

Gj,±:=^±^, (28) 

2Vr. a/v7 



are the Geljand-Dikii differential polynomials in u± := v^"^' it y^u^^', respectively. 
Proof. By identifying the coefficient of e^^+^ in (j22p we get 

Y^ JkU'-''^=0. (29) 

k+2l=2j+l 

From (I23p and by taking into account that {z — UcY = ^c + 4uc it is clear that only the terms 
JiU and J^U contribute to the coefficient of .^i in ((29]) . Thus we get the recursion 

relations 



(30) 

a,/3f ) = (^;,4 + 4ll(^)cfe+24^))/3fJr'^^ +4^c(2n(2)a^ + 4'Vfr^ 

which lead at once to the well-known third-order differential equation for the Gelfand-Dikii differ- 
ential polynomials 

a^G,-± = (^e4 + 4tx±6W + 2u±,^) G,_i,±. (31) 

D 

From ([28]) we have that Gi^± = =b-^n± and by using ([3T]) we find 

2 

G2,± = ^—= {vcU±;xx + 3u|) , 



2 

G3,± = ±^= (Wc^±,irx3!x + 10WcU±,5Sj;U± + hv^U^^-^ + lOu^) , 



yVc ^ 

+7{)VcUj^u±^xx + 70fcU±ii-|- X + 35n_|_ 



3 Strong conditions for critical points and the dou- 
bling property 

In what follows the notation F'^^' will represent the value of a function F at a critical point 
{xc,tc,Uc,Vc) of Q. We will also suppose that Vc 7^ 0. The following sets of critical points were 
considered in the applications of the Herniitian matrix model to quantum gravity |12j-|lUj. 

Definition 1. Given n >2 we denote by Cn the set of critical points of the hodograph system ([9]) 
which satisfy the ('strongj conditions 



{d^f)('^^ = {dudif)^'^=0, l<k<n, l<l<n-l, 
and such that ((a^+V)^'\ (^n^? /)^^^) / (0,0). 



(32) 



Due to ([To]) it is clear that C„ is also determined by the condition that all derivatives {d^ d\, fY'^' 
with (k, I) 7^ (1, 0) up to order n vanish and such that at least one n + 1-order derivative is different 
from zero. 

An alternative characterization of C„ can be formulated in terms of the integrals 



k 

'7 Vc 



Vz{z,t,) 



Ikitc,Uc,Vc) '■= ^ f '^^ "''fc+1/2' ' Jfc(*c, Uc, Vc) := ^. d) dz{z- Uc) ' ^fc"'^/^' • (^^) 



Indeed, Ik and Jk are proportional to {d^^^ f)^'^' and {dud^ fp'^' respectively. Hence it follows 
that 



Lemma 1. {xc,tc^Uc,Vc) G C„ if and only if 

h = 0, Jo = -2xc, 4 = Jfc = 0, forl<k<n-l , (/„, J„) / (0,0). 



(34) 



Let us consider the system of string equations ([5]) at points (x, tc, u, v) near a given critical 
point {xc,tc,Uc,Vc) G Cn- In terms of symmetric variables it reads 



dz 

2TTiz 



V,{z,tc)Y{z) =0, 



(35) 



VJ7 



dz 

27ri 



y,(z,te)U(z) 



-2x. 



Then if we set 



x = xe + e2"x, e = e2"+\ 



and assume (|25p we obtain a recursive method for determining the coefficients of w and v. Indeed 
([35]) is equivalent to 



^^ V,iz,tc)V^'^\z) =0, 



2TTiz 



f -^ V,iz, te) U^^'^ (z) = -2 xc Sjo -25j, 



(36) 



10 



For < J < n — 1 these equations are identically satisfied because of (|34p . For j = n we get from 
(f26]) - (|27|l that the equations (i36|) reduce to 

(37) 

or, equivalently, in terms of the Gel'fand-Dikii polynomials ()28p we obtain a pair of decoupled 
ordinary differential equations for u^ := v^"^' it y^u'^^ 

(jn±2^,lAGn,±iu±) = T^x, (38) 



Thus if the condition 

is satisfied, we have the so-called doubling property arising in the one-cut case of the Hermitian 
matrix model |10j-|llj. 

The first few cases of (1381) are 



\juvv i v^c Jvvv) \Vc u±7^ + o U-^j — 12 X, n — 2 ; 

{fLmv±^/v^fvuuv)\v'^u±^xxxx =120^, n = 3; 

\Juvvvv 31 \J V c J uvu-u-u j \V ^ U-i-^xxxxxx ~r -1-4 f^ U-|- U-i-^xxxx ~r ^oV^ U-^xxxU^^x ~r -^-L "^c ^5: 



za; 



+70 7;c -Ui Wi.xx + 70 7;c u± «! 5. + 35 u^ ] = 1680 x, n = 4. 
Finally for j > n the system (I36p yields the pair of equations 

± ( J. 7^^-^ + /. .f ^■)) = 0, t {J^ -?'' + ^' A^'^'O = 0' (39) 

which determine each pair (t(;(^(-'~"+^),t;(^(-'~"'+^)) recursively. 

4 Further critical points and regularized expansions 



Let us go back to the system (|TT|) for critical points of the hodograph equations ([9]) . It is equivalent 
to 

fv =0, fu = 2x, fuv = cr,/vfyy, 0" = ±1. (40) 

Let us consider solutions of the hodograph equations near critical points and assume that 

fy(tc,U,Vc) ^ 

(similar results are obtained by interchanging the roles of u and v). Then there exists an integer 
/ > 1 verifying 

{dt /,)(^) = 0, (0 < A; < /), {dl f,t^ + 0. (41) 

11 



As a consequence the first hodograph equation fv{tc, u,v) = can be used to eliminate the variable 
n as a function of v near {uc,Vc)- Indeed, from Weierstrass' preparation theorem it follows that 
near (uc, Vc) there exists a factorization 



f^{tc,u,v) = {Ao{tc,v) + ... + Ai_i{tc,v)u' ' + u')g{tc,u,v), (42) 

where Aj {0 < j < I) are analytic functions of v, and g is an analytic function of (n, v) which 
does not vanish near {uc,Vc)- Hence the solutions of the hodograph equation fy(tc,u,v) = near 
{uc, Vc) are given by the roots of the polynomial factor in (j42p 



Ao(*c, v) + ... + A-i(*c, v) n'-i + n^ = 0, (43) 

and consequently they are characterized by a Puiseux series 

u{tc,v) =Uc + 2_^ ak(tc){v -Vc)"^, (44) 

k>l 

for a certain integer 1 < m < I. If we now introduce the function 

H{tc,w):=fu{tc,u{tc,Vc + w"'),Vc + w"'), {w:={v-Vc)^) (45) 

then at t = tc the second hodograph equation in ([9|) reads H{tc, w) = 2x. 

Furthermore it is easy to see that as a consequence of the system (j40p we have 



On the other hand by differentiating the identity /„(tc, u(tc, Vc + w^),Vc + w"^) = we get 

Hence, we deduce that 

F(f) = 0. (46) 

In this way if we assume that there exists an integer n > 2 such that 

(dt H)^'^ = 0, (1 < A; < n), {dl F) « / 0, (47) 

then the hodograph equation H{tc,w) = 2x determines w as a function of x with a branch point 
of order n — 1 at x = Xc 

Definition 2. We will denote by Cimn-, (^ > 1, 1 < wi < /, n > 2) the set of critical points of ([9]) 
characterized by (|40|) -(|41| ) and such that 

1. u = u{tc, v) has branching order m — 1 at v = Vq- 

2. The corresponding function H{tc,w) defined by (|45p satisfies (j47p . 



12 



Example 



If we set t„ = for n 7^ 1, 3 then the system (jl0|) reads 

tci + 3tc3{ul + 2vc) = 0, Qtc3UcVc + Xc = 0, u^ = a ^/v^, a = ±1. (48) 

Let us consider the critical points {xc,tci,tc3,Uc,Vc) with tcs 7^ 0. They are given by 



9tc3^ 



Ur. = a ^ V. 



where (xctci^^cs) are constrained by the equation 



3Xr = 2(Tt, 



cl 



tcl \l/2 



9tc3 



As we are assuming that fc 7^ 0, we have that td 7^ and Uc 7^ so that (j4ip is satisfied by Z = 1. 
Moreover the hodograph equation 



leads to 



U{tc,w) 



tci+3tc3{u^ + 2v) =0, 



2[w- — ]] , ■u;:=t;-fr., ?^r. 



(49) 



t 



cl 



7^0. 



2yy ' 9t,3 

Furthermore H{tc,'w) := —12tc3u{tc,'w) (w + Vc) verifies H^wiic-,^) / 0. Therefore, it follows that 

(Xc-,tcl-,tc3-,Uc,Vc) e C112. 

The following statements will be useful for the subsequent discussion. They are easily proved 
by differentiating ([^5]) and the identity fv{t,u{w),Vc + w"^) = 0. 

Lemma 2. Given {xc,tc,Uc,Vc) G Cimn with Vc ^ then 

• If I = 1 then fyy does not vanish and 

H^l = -c7v-'/H3f^'i)+Avcf^il-4a^cfit)- 



(50) 



Ifl>2 then fi^ = ft^ = /iS^ = and 



Hi-) 



(c) (c) ^ Ac) 



Vc Juvv {Uw ) -\- ZVc JvvvUw ~T Juvvi JOT Vfl — i 



^ 'Jc J uw \ I^W j ) 



(51) 



for m >2, 



,W 



where Uw satisfies 



Vcjvvv\^w I \ ^JuvvUw ~r Jvvv — U, JOT m — i, 



^ Jvvv "-«! "J; 



(52) 



for 777- > 2. 



As a consequence we deduce the following conditions for critical points with ri = 2 
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Proposition 1. Given {xc,tc,Uc,Vc) G Cim2 with Vc 7^ 0, then it follows that 

1. If {xc,tc,Uc,Vc) eCii2, then 

Sa/i^) + AavJ^i - A^jit / 0. (53) 

2. If {xc,tc,Uc,Vc) G C112 with I > 2 then 

/ii'-t'c/iSlVo. (54) 

3. If {xc,tc,Uc,Vc) G Cim2 with l,m>2 then 

fit^O, and /ii = 0. (55) 



Proof. For / = 1 we have that ?tt, = 1 so that (j53p is a consequence of (I50p . For I > 2 and m = 1 
by substituting the solution of (j52]) into ([5T]) we obtain the condition ([5^ . 



Finally suppose that I > 2, m> 2. Then according to ([52]) and (JSTI) we have that Uw fvvv = 0, 
and Hww = Vc{uw Y fuvv Therefore the conditions for -ffiny / are fvvv = and fuvv / which 



proves ([55|) . 

D 

Let us consider the system of string equations ([5]) at points (x, t, n, v) near a given critical point 
{xc,tc,Uc,Vc) S Qmn- In terms of symmetric variables it reads 

(56) 

r dz 

i> --I4(z,t)U(z) =-2x. 

One may introduce stretched variables not only for x but also for t. The most symmetrical choice 
is 

x = xc + e'"x; tj=tc,j + e^''tj, j > 1; e = e^''^\ (57) 

Notice that edt^ = ed^. so that (fT2|) . (1571) and 

tf = Uc + ^ e-' u^-'^ , v = -Uc + ^ e-'u^-'^ , (58) 

j>2 j>2 

are consistent with the Lax equations ([2]). 

We will next concentrate on the case n = 2 and will provide a recursive method for determining 
expansions of the form (j58p near critical points in Cim2- Thus we set 

X = Xc + e X, tj = tcj + e ij, j > 1, e = e ' . 

From ([56]) , equating the coefficients of order P one finds 
• For < j < 3: 

— T4(z,t,)^ = 0, j—V,{z,QU^'^ = -2x,6,o, (59) 
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For j > 4: 



—-y,{z, u) — + i —-y,{z, t) 

2m z L, Zin z 



0, 



(60) 



Let us first analyze the system (|59|) . For j = it reduces to the hodograph equations 

Jv = U, /„ = Z Xc, 

and for j = 1 is triviahy verified since U = V =0. Furthermore, it is straightforward to see 
that for J = 2 both equations in ([59]) reduce to 



(61) 



Hence to proceed further it is required to distinguish the two types of critical points corresponding 
to Cii2 and Cim2 with I > 2. 

Case Cii2 

In this case fvv / and ()6111 implies 

_(2) _(2) ^ -(2) 



For i = 3 both equations in (I59p lead to 



^ZJ(3)=_^^(3). 



By setting j = 4 in (j60p we find that 



tf 



<4) = _^(^(4) + l^(2)+ 1 -(2)^ 



a fy{t,Uc,Vc 

\/'Vc Jvv 



+ 



'J Jvv 



Juvv ^ V ^c /, 



vvv 



vU + -I^(^)^ 



Moreover v^"^' = —a^/v^u^'^' must satisfy the Painleve I equation 

r^'^ (^2 + -^^'^') = 6 (2 X - fu{t, Uc, Vc) + aV^cfvil u„ '. 
where 



rW := -a ^(3 /i^) + ^vc f^t - 4a V^/l^). 
Notice that since (xc, to^o'^c) G C112 then Proposition 1 implies F^'^^ 7^ 0. 



(62) 
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To proceed further one uses induction in the recurrence system for U , V to prove that the 
only coefficients in the expansions of C/ , V depending on v^^', v^^~^\ u'-'~^\ w^^', w^^~^' and 
w^^""^' are those corresponding to / = 1 or 2. Moreover one finds 

r?p') = 2l;(^-) + ^I;(4-')+2lZj(2)uj(^'-2)+Di^'\ 

where A\ , B^ , C^ and D^ , j > 5,i = 1,2 are differential polynomials in I;(^\ . . . ,v^^^'^\w^^' , . . . ,w^^^^\ 
Then, for j > 5 substituting (|24p into the first equation in (j60p yields expressions 



where -ftTj are differential polynomials in v^"^' , . . . ,v^^ ^'. Moreover, both equations in (j60p imply 
that tJ"~^^ must verify a second order linear differential equation of the form 

r^'^) (tJ^r^) + l^u(2)^0-2)) = i7^._2(t,^;(2), . . . ,^0--3))^ (63) 

where -ffj-2 are differential polynomials in v^'^\...,v^^~'^K In this way, we have a scheme for deter- 
mining the coefficients u"-*, w^^' in (J14p . 

We notice that if {v^'"' , / > 2} is a solution of ([62]) - (f63l) {j > 5), then as a consequence of 
the symmetry transformation (J16p - (|19p we have that {(— 1)%' ', / > 2} is also a solution of these 
equations. Consequently, the differential polynomials H21+1, / > 1, in the right hand side of (|63p 
are odd polynomials in v^'^^'^^', 1 < i < / — 1 and their derivatives (for example H^it,^'^') = 0). 
Hence we may set iJv^J+i) = for all j > 1. Analogously, since the equations ([62|) - ([63P can be 
written in terms of {w^', I > 2}, we may set w^'^^^^' = for all j > 1. 
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Case Cim2, I > 2 

In this case the system (j59p is trivially satisfied. Moreover, for j = 4 we get from (j60p that 

'A^') {vJJ^^ + 3I;(2)' + 37;,IZJ(2)') = 6 (/(^ (/„(t,ne,t^c) - 2x) - t;e/iSl A(t,u„T;,)), 

(64) 
, A(^) (t;cIIjS + 6lZj(2)u(2)) = 6 (/(^ /,(t, ue, t;e) - f^t Unit, n„ ^,) - 2 x)) , 

where we are denoting 

From Proposition 1 we have that A*-"^) ^ for critical points in Cijn2 with I > 2. 

In terms of the variables u^ ' = u'^' it y^u'^-*, the system ()64p decouples into the two Painleve 
I equations 

ifivv ± \/^^/iSl)('^c U±;^ + 3 (n±)^) = 6 (2 X - fu{t, Uc, Vc) =F V^fv{t, Uc, Vc)). 

Thus the doubling property is satisfied in this case. 

In general, for j > 5 (j60p leads to a second order linear system for w^^~'^\ v^^~'^> of the form: 

2v, (Ijgr'^ + 6 w;(2)w;(.'--2)) + l2u(2)i;{i-2) ^ M,_2(t,IZj(2), . . . ,i;(^--3),^(2)^ . . . ^^0-3))^ 

(65) 
vMi^^'^ + 6lZj(2)uO--2) + 6u(2)lz;(^-^2) ^ Ar^._2(t,iu(2), . . . ,iuO-3),u(2), . . . ^^O'-s)). 

where Mj_2 and A''j_2 are differential polynomials in w^"^', . . . ,w^^^^' and m \ . . . ^v^^K Thus, 
we have a recursive procedure to construct the coefficients IZJ"\ v^^', j > 2. 

Again, due to the symmetry ([TUD-([THD, it is clear that if {w^^', v^^'), (j > 2) is solution of ([651) 
then ((— l)-'uJ"\ {—lyv^^'), {j > 2) is also a solution of (j65p . Hence for odd j the differential 
polynomials Mj-2, Nj-2 are odd polynomials in Itj'2'+^), u^^'+ii^ I > i and their x-derivatives. 
Therefore, we may set l(;(2'+i) = Q, if(2'+i) = Q for all / > 1. 

4.1 Critical processes in ideal Hele-Shaw flows 

In view of the properties of asymptotic solutions of the KdV equation [I3]-[T5] and the NLS equation 
[IB] it should be expected that the inner expansions provided by the double scaling method will 
be relevant when the solutions of the dispersionless or the dispersionful Toda hierarchies reach a 
point of gradient catastrophe. We next consider an application to an ideal model of Hele-Shaw 
flows supplied by the Toda hierarchy. 

A Hele-Shaw cell is a narrow gap between two plates filled with two fluids: say oil surrounding 
one or several bubbles of air. In the set-up considered in |18) (see also |19j-|20]) air is injected in two 
flxed points of a simply-connected air bubble making the bubble break into two emergent bubbles. 
Before the break-off the interface oil-air remains free of cusp-like singularities and develops a smooth 
neck. The reversed evolution describes the merging of two bubbles. The analysis of [18] concludes 
that before the merging, the local structure of a small part of the interface containing the tips of 
the bubbles is described by a curve Y = Y{X) which falls into universal classes characterized by 
two even integers (4n, 2), n > 1, and a finite number 2n of deformation parameters t^. Assuming 
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symmetry of the curve with respect to the X-axis, the general solution for the curve in the (4 n, 2) 
class is 

y(z) := (^%£±llii^) V(--a)(z-6), X = z. (66) 

^ ^J [z — a){z — h) ^ffi 

where o and h are the positions of the bubbles tips. Due to the physical assumptions of the problem, 
the expansion 

2n oo -j^ 

Y{z) = Y,{k + l)tu+iz^ + Y,^, z^oo, (67) 

fc=l fc=0 

must satisfy the conditions Yq = t (physical time) and Yi = (bubble merging condition) which 
determine the positions a, h of the tips. However if Yi / the evolution process leads to a critical 
point in which cusp-like singularities appear. 

As it was shown in [H] the positions of the bubbles tips are determined by the pair of hodograph 
equations 

oo oo 

^fetfcrfc_i(u,z;) = 0, ^^ktkrkiu^v) + 2x = {), (68) 

fe=i fc=i 

where t^ = for A; > 2n, Yi = 2x and 

r:= ,, ", ==E^%^' a:=«-2V^, h:=u + 2V^. (69) 

These are precisely the hodograph equations ^. Thus the double scaling limit method can be used 
to regularize the solutions of ()66p -()68|) at critical points in terms of inner expansions of solutions 
of®. 

As an example let us analyze the critical process of a merging of two bubbles studied in section 
VII of [H] . We set t := ti, ^2 = 0, t^ = 0, n > 3 and fix t-^ to a given constant value c. Thus the 
Hele-Shaw interface is locally characterized by the curve 



Y(X) = 2,c(X + u)^ {X - uf - 4.V, (70) 

and ([68|) reduces to 

i + 3c(u^ + 2w) = 0, Qcvu + x = Q. (71) 

This is the hodograph system (|48p and its solution is given by 



t2 



22/3 9 c \/9 V81 c2 ^4 + 4 (^^273 _ 81 c a;2 _ 2 13 18 c 



\/9 \/8rc2 x^ + 4ct3:E"2 - 81 cx2 - 2t3 

H — TT^ , ■'^ 



(72) 



18v^c ' Qcv' 



As we have seen above the critical points {xc-,tc,tci,Uc,Vc^ with tc 7^ are in C112 so that our 
scheme can be applied to provide a dispersive regularization of (f72|) near critical points. 

Setting X = Xc in ([72]) one finds that the outer approximations for u and v near t = tc are given 
by 



w ~ ^io«t(i) :='«c - oV -(*c-i), w ~ ^^o«t(t) := ^c + YW-(ic-i), as t ^ t^ , (73) 



where Uc = cr ^Jvl. Since we are considering a critical point in C112 we introduce the stretched 
variables 



and consider the inner expansions for u and v 

u '-^ Uin := Uc + e u^ ' {x — Uct) , v ^^ Vin{t) := Vc + e v^ '{x — Uct), as t —i- 0, (74) 
where w'^' = —v^'^'/uc and v^'^' verifies the Painleve I equation 

Ilg + ^(l7(^))^ = ^(x-..t). (75) 

The inner approximation at a; = must match the outer approximation in an overlap interval 
which has both t — tc small and t large. Writing the outer approximations (|73p in terms of the 
inner variable t we have 

Uout{t) = Uc-e^ -\ — , Vout{t) =Vc + e^ -^ 
3 V c 3 

Hence, it is clear that matching requires a solution w"^' (s — Uc t) of (|75p satisfying 

t;(^)(-nct) ~ ^A/ — , as t ^ -00. (76) 

O V c 

Now, if we introduce the change of variables 
it follows that W must satisfy the P-I equation 




and 



y/ n^ —\ -^ as ^^00, 
V 6 



so that W must be the tritronquee solution of the P-I equation [2l]-[25]- Using a numerical 
approximation of W one finds [27] that the regularized evolution is as follows: the right bubble 
develops a cusp, then a new bubble appears at this cusp and it grows until it merges with the tip 
of the left bubble. Finally, the tip of the left bubble absorbs the new one and stays joined to the 
right bubble. 
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